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gauge-invariant matrix models. Special attention is paid to solving the models at 
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1 Introduction 



Matrix models are usually associated ^ with discretized random surfaces (or string the- 
ory) in ad < 1-dimensional embedding space and, in particular, with 2D quantum gravity. 
The simplest Hermitean one-matrix model corresponds p|, Q to pure 2D gravity while a 
chain of Hermitean matrices describes |Q, ^ ^ 2D gravity interacting with d < 1 mat- 
ter. A natural multi-dimensional extension of this construction is associated with 
induced lattice gauge theories. The matrix models which describe induced QCD can be 
constructed in a similar way 0. 



1.1 Random surfaces 

The typical problems which reduce to matrix models are associated with a statistical 
ensemble of random surfaces whose partition function is defined generically as 

Zrs = T.^-''^^'^- (1-1) 

s 

Here A{S) is the area of the surface S which can be either closed or open and a stands 
for the string tension. 

There is a lot of examples of such systems in quantum field theory: 

• Strings either fundamental (graviton) or secondary (hadrons). 

• 3D Ising model (the boundary between different phases is two-dimensional). 

• Lattice gauge theory at strong coupling. 

• 1 /A^-expansion of QCD §. 

• 2D quantum gravity. 

The last system is described by the Euclidean partition function 

Z^o = J Dge-I'^'^'M^-^'^) = J Dg I '"''^^''+0 (1.2) 
where A stands for the cosmological constant and x is the Euler characteristics of the 2D 



world. The path integral in Eq. (1.2) is over all metrics gnu{x 



1.2 Dynamical triangulation 

The idea of dynamical triangulation of random surfaces is to approximate the surface by 
a set of equilateral triangles. The coordination number (the number of triangles meeting 
at a vertex) is not necessarily equal to six which is associated with internal curvature of 
the surface. The partition function (|1.2| ) is approximated by |jll| 

^DT = Ee^^'"'^Ee-^"* (1.3) 

where one splits the sum over all the triangles into the sum over the genera, g, and the 
sum over all possible triangulations, T^, at fixed genus g. Remember that (? = for a 
sphere, g = 1 for a torus and x = 2(1 — g). In ( |1.3| ) rit stands for the number of triangles 
which is not fixed and is a dynamical variable. 

The exponential suppression with nt provides the convergence of the sum over Tg 
in (|1.3| ) at least for large enough A. However, the sum can diverge for some values of A 
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Figure 1: A graph constructed from equilateral triangles (depicted in the bold lines) and 
its dual one (depicted in the thin lines). 

due to the entropy factor (the number of graphs). It is crucial for what follows that the 
total number of graphs of genus g with n triangles grows at large n as a power of n [|10| : 

^ 5{nt ~n)= e^^'^n-^^ (l + 0{n~^)) , (1.4) 

where Ac does not depend on g. For this reason the genus g contributions to the string 
susceptibility 

/ = ^,ZnT - E e ^(^-^)(A - Ac)-^« , 7. = -fe. + 3 (1.5) 

simultaneously diverge as A — > Ac + 0. This is the point where the continuum limit is 
reached and the discrete partition function Zdt approaches the continuum one Z2d- 

A similar dynamical triangulation can be written [|l| for the partition function (|1 . 1|) 
of 2D surfaces embedded in a (i-dimensional space. The case of 2D gravity is associated 
with (i = 0. 

1.3 Large-A^ matrix models 

The partition function ( |1.3| ) can be represented as a matrix model. The dual graph for 
the set of triangles coincides with the graph in a = quantum field theory with a 
cubic interaction as is depicted in Fig. |T]. The precise statement is that Z^t equals to the 
partition function of the N x N Hermitean one-matrix model 

Z,,,^ e^'^ = |t/$e-^*^^(*) (1.6) 

with N = exp(l/G) and the cubic coupling constant ^3 = exp(— Ac). The integration 
measure in Eq. (|1.6|) is 

N N 

= '[[dRe^ijdlm^ijY[d^ii (1.7) 

i>j 1=1 

and V{^) = + ^3$^ is a cubic potential. 
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The fact that Z/jt = Zim can be proven analyzing the "fat-graph" expansion of (|1.6|) 
with the propagator 



which leads for log Zim to the factor A^^ associated 0] with a graph of genus g. 
The partition function ( p..6| ) with the general potential 

oo 

V{^) = J2t,^^ (1.9) 

j=0 

is associated with a discretization by regular polygons with j > 3 vertices whose area is 
j-2 times the area of the equilateral triangle. 

1.4 Double-scaling limit 



A question arises how the system described by the partition function ( |1.6|) can undergo 



a phase transition for A ^ which is associated, as is discussed in Subsect. |1.2| , with 



the continuum limit. While the system is at d = 0, a (third-order) phase transition of the 



Gross-Witten type [|1T| is possible as A^ oo when the number of degrees of freedom 
becomes infinite. Therefore, the continuum limit is reached at A^ ^ oo and A Ac. 

The N = oo limit corresponds to planar diagrams or genus zero (the spherical approx- 
imation). Higher genera are suppressed as N~^^. 

One can utilize, however, the fact that •jg, which is defined by Eq. ( [L.5|) , linearly 



depends on g [12 



7. = 2 + ^(^?-l). (1.10) 



Therefore, the parameter of the genus expansion near the critical point is 

1 
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Ar2(A- Ac)5 



;i.ii) 



and can be made finite if (A — Ac) ~ N's as N ^ oo 0. This special limit when the 
couplings reach critical values in a A^-dependent way as A^ ^ oo is called the double 
scaling limit. The double scaling limit of the Hermitean one-matrix model allowed |Q to 
construct the genus expansion of 2D quantum gravity. 

1.5 The Kazakov— Migdal model and induced QCD 

A natural d > 1-dimensional extension of ( |1.6| ) is the Kazakov-Migdal model which is 
defined by the partition function 

Zkm= [[dU^{x)[[d<^^e "V >^~^ J. (1.12) 

X,fJ, X 

Here the integration over the gauge field ?7^(x) is over the Haar measure on SU{N) at 
each link of a (i-dimensional lattice with x labeling its sites. The model ( |1.12| ) obviously 
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recovers the standard d < 1 matrix chain if the lattice is just a one-dimensional sequence 
of points for which the gauge field can be absorbed by a unitary transformation of 

The large-iV solution of the Kazakov-Migdal model in the strong coupling phase 
|l^-|T^ is associated [0 with an unbroken extra Z{N) symmetry of the partition func- 



tion ( |1.12[ ) and, therefore, with infinite string tension (see []TB[ for a review). One can 
easily modify the model in order to have a phase transition (with decreasing the bare mass 
parameter) after which the Z{N) symmetry is broken in some sense and normal area law 
associated with finite string tension is restored. While the arguments are based on the 
mean field analysis, they look quite reasonable because the phase transition occurs only 
for systems which are not asymptotically free. 

The simplest model of this type is the adjoint fermion model which is defined by the 
partition function M 



'AFM 



I nrff/,(^)n^^^^^x.e-^-[™] (1.13) 



where '^x and are the N x N matrices whose elements are independent anticommuting 
Grassmann variables, Sf[^,^,U] is the lattice fermion action 



Sf = J2NtT {Vpi^x^x) -J2[^xP-U,{x)^x+,Ul{x) + ^x+,P^Ul{x)^xU,{x)]) (1.14) 



D 

I 

11=1 

and are the standard projectors. 



2 Higher genera in one-matrix model 

The Hermitean one-matrix model was first solved in genus zero in Ref. [|T^ by the method 
of the saddle-point integral equations for the spectral density. The more powerful orthog- 
onal polynomial technique allowed to calculate the partition function up to genus two for 
the quartic interaction I review in this section the method of solving the Hermitean 
one-matrix model which is based on the loop equations and provides an algorithm for 



genus by genus calculations. The explicit results El], on calculation of the partition 



function with an arbitrary potential and all correlators are presented up to genus two. 
2.1 Loop equation 



All correlators of the Hermitean one-matrix model ( |1.6|) can be obtained from the Laplace 
image of the Wilson loop 

WiX) = (%T^) (2.1) 



A-<l> 

where the averaging is w.r.t. the same measure as in (|1.6| ). As is explained in Ref. 



W{X) is associated with the sum over discretized open surfaces with one fixed boundary. 

The correlator (|2.1| ) can be obtained from the free energy, F, by applying the loop 
insertion operator, d/dV{X): 
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W{\) is determined by the loop equation (see for a review 



, P:P^Wiu:) = W\X) + ^^-^WiX) (2.3) 

where the contour Ci encloses counterclockwise singularities of W{u). The contour inte- 
gration acts as a projector picking up negative powers of A. The second term on the r.h.s. 
of the loop equation ( p.3|) is expressed via W{X), so that Eq. (|2.3| ) is closed and unam- 
biguously determines W{X) imposing the boundary condition XW{X) ^ 1 as A — > oo. 
The genus expansion of W{X) and of the free energy, F, is defined by 

mX)-Zj^W,W, F^^-F. with MMA,^-^. (2.4) 
2.2 Genus zero solution 

To leading order in l/N"^ one can disregard the second term on the r.h.s. of Eq. ( |2.3| ) 
which reduces to the quadratic equation 

^rMw(.) = WHX). (2,5) 



The one-cut solution to Eq. ( p.5|) reads p3| 

du V'{u) 



Wo{X) 



{X-x)iX-y) 



Ci 4:7ii X — uj\^ (lu — x)(uj — y) 



(2.6) 



where x and y are determined by 

duj V'{iu) ^ f du iuV'iiu) 



2712 ^(cu-x)(cu-i/) ' -^^1 2vrz ^(cj - x)(cj - y) 



(2.7) 



The formulas (|2.6|) , ( p.7|) solves Eq. ( p.5|) when V'{uj) is polynomial (which is usually 
associated with the one- matrix model) or has singularities outside the cut [y, x] . 

Doing the contour integral in (p.6|) by taking the residues at = A and u = oo, one 
finds 

WoiX) = \ \V'{X) - M{X)^{X-x){X-y)] (2.8) 

where M(A) is a polynomial in A of degree J-2 if ^(A) is that of degree J. This form of the 
genus zero solution is convenient to determine the spectral density, p(A), which describes 
the distribution of eigenvalues of the matrix $ at the large- saddle point. Calculating 
the discontinuity of W{X) across the cut [?/, x], one gets 

p(A) = Im W{X) = -M{X)^{X-y){x-X) X G [y, x] . (2.9) 

Eq. ( p.7|) guarantees that / dX p(A) = 1. 

The spectral density p(A) given by ( p. 91) vanishes under normal circumstances as a 
square root at both ends of its support. The critical behavior emerges when some of the 
roots of Af (A) approach the end point x or y. 
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2.3 The iterative procedure 



The iterative procedure of solving the loop equation is based on the genus zero solu- 
tion ([^.61). Inserting the genus expansion ( |2.4| ) in Eq. (|2.3|), one gets the following equation 
for Wg{\) at c/ > 1: 

L ^ r^^^^^^ - 21^o(A) W,{X) = Y.W.W W.-.'W + , (2.10) 

which expresses Wg{X) entirely in terms of Wg/{X) with g' < g. This makes it possible to 
solve Eq. ( |2.10| ) iteratively genus by genus. 

The iterative procedure simplifies if one introduces, instead of the coupling constants 
tj, the moments Mk and Jk defined for /c > 1 by 

/ V'H I - f ^'^^^ (211) 

These moments depend on the coupling constants t/s both explicitly and via x and y 
which are determined by Eq. (|2.7| ). Notice that Mk and Jk depend explicitly only on tj 
with i >k + 1. 

The main motivation for introducing the moments ( p.ll| ) is that Wg{X) depends only 



on 2 X (3(7 — 1) lower moments ( 2 x (3(7 — 2) for Fg ) This is in contrast to the 

t-dependence of Wg and Fg which always depend on the infinite set of tj's (1 < j < oo). 



2.4 Genus one and two results 



To find Fg, one first solves Eq. ( |2.10|) for Wg{X) and then uses the last equation in 
The result in genus one reads pl| 

— hid 
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1 1 

Fi = InMi In Ji 

24 24 



(2.12) 



where d = x — y. 

An analogous calculation in genus two yields p2 



119 
7680 Ji^d* 



119 
7680Afi^d4 



181 J2 

480 Ji^d^ 



181 M2 



480 Ml 



3 J2 

64Ji^Mi(i3 



3M2 
64 Ji Mi^ d? 



IIM2Z 



'40 Ji^d^ 40Mi*d2 192Mi^d2 



43Af3 



+ 



+ 



2IJ23 _ 29 J2 J3 



35 J4 



WOJi^d 128 Ji^ d ~'~ 384 Ji^d 160Mi^d ' 128 Mi* d 384Mi^(i 

Some of these coefficients have an interpretation in terms of the intersection indices on 
moduli space |2^ , the others are associated with characteristics of the discretized moduli 
space . 

Since Fg is known, the genus g contribution to any connected correlator 



43 J3 
192 Ji^d2 

21 M2^ 



+ 



J2 M2 



64JP~mP^ 



17 

128 Ji Ml d* 



+ 



29 M2 M3 _ 35 M4 



(2.13) 



tr 
N 



tr 



2-2s 



d 



d 



dti, ' ' ' dti 



(2.14) 



can be calculated by the differentiation. It depends on at most 2 x {3g — 2 + s) lower 
moments []2ll , p2| . 
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To obtain explicit formulas, say, for the symmetric quartic potential when all tj = 
except t2 and t^, one should solve Eq. ( ^^7]) for x = —y: 



X 



/('M +A 

3t4 V V3t4/ 3t4 



(2.15) 



and express the moments ( |2.11D via t2 and which is given by algebraic formulas 



Ml = Ji = 2t2 + 6Ux^ , M2 = -J2 = 8Ux , M^ = 



At 
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Mk = Jk = for A; > 4 . 

(2.16) 

The results for Fi and F2 in the case of the quartic potential are in agreement with those 



ofRef. 20 



3 Solving matrix models at d > 1 

The Kazakov-Migdal model was originally studied at large- by the Riemann-Hilbert 
method [1^. The explicit solution for the quadratic potential [1^] was reproduced [1^ 
by the loop equations. The equivalence of the two methods was shown for an arbitrary 



potential in Ref. [|23 where the relation to the Hermitean two-matrix model was utilized. 
This approach allowed to solve explicitly the Kazakov-Migdal model with a logarithmic 
potential and the adjoint fermion model with the quadratic potential p9[ in the strong 
coupling phase. 



3.1 Loop equation for one-link correlator 



Let us define for the Kazakov-Migdal model (|1.12|) the loop average and the one-link 
correlator, respectively, by 



W{\) 



tr 



(3.1) 



The definition of W{X) is similar to Eq. ( [^.ID while (?(//, A), which is symmetric in u and 
A due to invariance of the Haar measure, dU, under the transformation U W , is absent 
in the one-matrix model. Expanding (j(z/. A) in 1/z/, one gets 



n=l 



tr 
N 



{^:u,{x)^-L-uli 



(3.2) 



The correlator G{h', A) obeys in the large- limit the following equation ^1 

[ —^^^G{uj,X) = W{u)G{u,X) + XG{u,X)-W{iy), (3.3) 

JCi 2711 U — UJ 

where the contour Gi encircles counterclockwise the cut (or cuts) of the function G(a;, A) 
and 

V\lo) = V'{uj) - {2d - l)F(cj) . (3.4) 

The function 



tr 



Fiuj) = Y.F^u\ Fo= j^i^-EFn^ 



n=0 



n=l 



(3.5) 
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is determined by the pair correlator of the gauge fields 



tr 



^ - (r<|.") 

n=l 



(3.6) 



where $ and \1/ play the role of external fields and {a = 1, . . . , iV^-1) stand for the 
generators of the SU{N). Eq. (|3.6|) holds [jl3|, [l^ at = oo. The choice of Fq which is 
not determined by Eq. (p.6|) is a matter of convenience p8|. 



Taking the 1/A term of the expansion of Eq. (|0| ) in A and using Eqs. ( p.4|) and ( |3.6|) , 
one arrives at the equation for iy(z/) which coincides with Eq. ( |2.5D for the Hermitean 
one-matrix model where V is substituted by 



r{X) = V'(A) - F(A) . 



(3.7) 



The potential V{uj) is, generally speaking, non-polynomial and has singularities on the 
complex plane outside of the cut (or cuts) of W{uj). 

It is worth mentioning that Eq. (|3.3|) coincides with the loop equation for the Her- 



mitean two- matrix model |30]. This is because at d = 1/2, which is associated with the 
Hermitean two-matrix model, the last term on the r.h.s. of Eq. ( |3.4D disappears and one 
gets just V{uj) = V{u). 



3.2 The master field equation 



To analyze the model (|1.12| ), let us consider the Hermitean two-matrix model with the 
potential 



v($) = E 



m=l 



m 



The solution for W{X) versus V(A) is determined by the equation |27 

E 9mGm-l{X) = XW{X) - 1 



(3.8) 



(3.9) 



m>l 



which is just the l/u term of the expansion of Eq. ( p.3|) in l/u. 

The functions G'„(A) are expressed via W{X) using the recurrence relation 



Gn+l{X) 



duo V{uj) 
Ci 27ri X — uo 



Gn{^) - WiX) Gn{X) , Go(A) = WiX) 



(3.10) 



which is obtained expanding Eq. (|3.3|) in 1/A. If V(A) is a polynomial of degree J, Eq. (|3.9|) 
contains W{X) up to degree J and the solution is algebraic pO| . 



As is proven in Ref. p7 



]1) are 



i) Equations which appear from the next terms of the 1 / z^-expansion of Eq. 
automatically satisfied as a consequence of Eqs. (|3.9|) and ( p.lO| ). 

ii) (j(z/. A) is symmetric in and A for any solution of Eq. (|3.9|) . The symmetry re- 
quirement can be used directly to determine W{X) alternatively to Eq. ( p.9|) . 

Since the approach based on Eq. ( ^.3[ ) is equivalent |^ to that of Ref. |T^, G(z/, A) 
can be expressed via W{X) as follows [|1^, [1^, |31| 



G{u, X) = l-expUl^ 



dio 



Ci 27ri u 



UJ 



log (A - r±(^)) 



(3.11) 
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where 



?^±(A) = ±t7rp{X) 



F{\) + W{\) 



(3.12) 



(3.13) 



This relation in ci = 1 was advocated in Ref. ||3^ studying the large- asymptotics of the 
integral over the unitary group in ( |1.12D (the Itzykson-Zuber integral). Since G(z/, A) is 



The condition for the r.h.s. of Eq. ( p.ll| ) to be symmetric in u and A is 

r±(r^(A)) = A. 



symmetric in u and A, the master field equation |]T3 

duj 



W{\) = ±( -^log(A-r±(c.)) 



(3.14) 



obtained as the 1/v term of Eq. ( p.ll|) , will be satisfied as a consequence of Eq. ( |3.13|) 
which guarantees the symmetry. 

All three equations (p.9|), ( p.l3| ) and ( p. 141 ) seem to be equivalent. It is a matter of 
practical convenience which equation to solve. The only known explicit solutions exist for 
the quadratic potential |TB[ and the logarithmic potential . 



3.3 Explicit solution for logarithmic potential 

Let us choose the following potential of the Hermitean two-matrix model 



V($) 



-{ah + c) log (6 - $) - a<l> = ^ — + t$ 



m=2 



(3.15) 



where a and h are real. In the one-matrix case this potential is associated with the Penner 
model 



The quadratic potential is recovered in the limit 
a 



a,b ^ oo , 
The solution to Eq. 



c ~ 1 



(quadratic potential) . 



(3.16) 



for G{i', A) versus W{X) with V{uj) given by Eq. ( p.l5| ) is 



{a+X)W{X)-l 
b-u 



(3.17) 



W{X) is determined by Eqs. (|3.9|), ( p.lO| ) which reduce to the quadratic equation for W^u) 
of the form of Eq. (|2.5| ) for the Hermitean one-matrix model with the logarithmic potential 

\7($) = -(ab + c) log (6 - $) + (a6 + c + 1) log (a + $) - (a + 6)$ . (3.18) 



G(z/, A) given by Eq. (|3.17|) is indeed symmetric in u and A providing Eq. ( p.5|) with V = V 
is satisfied. 

Since V'{X) is known, the function -F(A) can be determined from Eq. ( |3.7| ) to be 



F(A) 



6A - c - 1 
a + X 



(3.19) 
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Now Eq. ( p.4|) determines the potential 

y(<|.) = -(a6 + c)log(6-<l>)-(2c/-l)(a6 + c+l)log(a + <l>) + [(2c/-l)6-a]$. (3.20) 

This formula recovers at d = 1/2 the potential ( p.l5| ) of the two- matrix model and at 
d = the potential ( ^.18 ) of the associated one-matrix model. 



In the Gaussian limit ( |3.16| ) all the formulas recover the ones for the solution found 
by Gross |T^. For this reason the one-cut solution (|2.6| ), (|2.7|) with the potential ( p.l8|) is 



always realized for iy(A) if a and b are large enough for the points b and —a to lie outside 
of the cut. 

While the above solution for the logarithmic potential was originally obtained 
solving Eq. (|3.3| ), it is easy to show that it satisfies Eq. (p.l3|) and, therefore, the master 
field equation ( |3.14| ). Let us notice for this purpose that r±{X) given by Eq. ( |3.12| ) with 
our solution for V'(A), F{X) and W{\) satisfy the following equation 

P(r±(A),A)=P(A,r^(A)) = (3.21) 

with 

V{u, A) = -X^u^ + {b- a)Xu{X + u) + ab{X^ + v^) - Xv[o? + 6^ + 20+1) 

+ (A + z/)(6-ac + 6c) -c2-62 + (a6 + c)[(a + 6)iy(6) -1], (3.22) 



where the constant Wi^) depends on the type of the solution of Eq. ( |2.5| ) (one-cut or 
more-than-one-cut solutions). Since this Div^X) is symmetric in v and A, Eq. ( |3.21|) 
implies that Eq. (|3.13|) is satisfied. 



According to Ref. [^], the pair correlator of \J and f/^ for the potential ( |3.20| ) can be 
calculated taking the discontinuity of ( |3.17|) both in v and A across the cut (cuts). For 



the solution ( p.l7| ) one gets 



C{v, A) = , \ Disc, Disc, G{v, A) = ^^^^""^^^^ , ^, A e cut (3.23) 

with I>(z/, A) given by Eq. (|3.22| ). In the Gaussian limit (|3.16|) when Wi})) 1/6, one 
recovers the result for C(z/, A) in the case of the quadratic potential. 

Since p(i/) is real for u 6 cut, the roots of the denominator, r±{i>), are complex so that 
C(z/, A) has no singularities at the cut. The arguments of Ref. ^4] suggest, therefore, that 
the Kazakov-Migdal model with the logarithmic potential ( p. 20 ) always remains in the 
phase with infinite string tension. 

The potential ( |3.20| ) of the d- dimensional matrix model admits the "naive" continuum 
limit when 

$ = e^-V , V{^) = e'^v{(j)) (3.24) 
with and f (0) being finite as £ — > 0. This continuum limit is reached providing 



as e^O (3.25) 



and results in a d-dimensional continuum quartic action |^ . The procedure of taking the 
"naive" continuum limit works for < 4 where b given by Eq. (|3.25|) is divergent. This 
is precisely where the scalar theory with the quartic interaction is renormalizable. Thus, 
one can look at the logarithmic potential ( p.20| ) as at a latticization of the quartic one for 
d<A. 
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3.4 Fermionic matrix models 



The correlators of arbitrary powers of '^x'^x at the same site x are determined for the 
fermionic matrix model ( |1.13| ) by 



tr 
N 



( 

\jj2 _ ^^^^ 



(3.26) 



d = 0) reads p 



The analogue of Eq. ( |2.3| ) for the fermionic one-matrix model (defined by ( |1.13| ) with 

-W{X). 



(3.27) 



This equation is identical to Eq. ( p.3|) for the Hermitean one-matrix model with the 
logarithmic potential 

V($) = 1/^($) +21og<l> (3.28) 

and $ = These two models are equivalent to all orders of the 1/iV-expansion. 

However, the genus expansion has now alternating signs and is convergent contrary to 
Ref. i. 

Let us define for the fermionic matrix model ( |1.13D on a d-dimensional lattice the 
odd-odd and even-even one-link correlators: 



W{u, X) 



7^ ( .^-'f.vi./ ^^^) 



A 



x+n ^ x+fi 



-"^x+.Uli 



^> ^ X + fl^ x+^ 



(3.29) 



where -l-(— ) are associated with the positive (negative) direction fi. 



The fermionic analogue of Eq. ( p.3|) consists of two equations p9 

^Y1Mg{uj,X) = Wo{iy)G{u,X)+\XW{iy,X)-Wo{iy) 
Ci Ani u — CO L 

^Ylit^W{uj,X) = Wo{u)W{u,X)-2-W{iy,X)-XG{u,X). 
Ci 4:7ii V — uj V 



(3.30) 



The solution to Eq. (|3.3CI|) for the quadratic potential Vp\ 



iS) = 17100"^ is 



Wo{X) = - 



/iA + ^-iV/x2A4 + 4 



with 



/i 



{D - l)m + DJm? + A{2D - 1) 



(3.31) 



(3.32) 



{2D - 1) 

This solution agrees ^ with the result ||3^ for lattice QCD with fundamental fermions 
at vanishing plaquette term. 
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4 Conclusions 



There is no problems to calculate higher orders of genus expansion in the Hermitean 
one-matrix model using the iterative method which is described in Sect. 2. It would be 
interesting to perform analogous calculations for the potential (p.28|) which is equivalent 



to the fermionic model whose genus expansion is expected to be convergent since the 
integral over the Grassmann variables in ( |1.13|) converges. 



The method of solving the Kazakov-Migdal model which is described in Sect. 3 reduces 
it at large- to the Hermitean two-matrix model. This goes along with Ref. ||3^ where 
the conformal field theories in < 1 are obtained from the two-matrix model. For d > 1 
the potential V of the two-matrix model should be, presumably, non-polynomial rather 
than polynomial as for d < 1. 

The approach of Sect. 3 works, however, only when singularities of V lie outside of 
the cut. This is not the case for the solution of the Riemann-Hilbert equations found 



by Migdal |T3[ which exhibits a non-trivial critical behavior. While this solution is not 
associated with induced QCD, it would be interesting to find out what physical system it 
corresponds to. 

The discussed solution of the d > 1 models are associated with the strong coupling 
phase, i.e. the phase with infinite string tension. An open question is whether the de- 
scribed approach can be extended to the phase with area law where QCD is induced. 
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